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1.  INTRODUCTION 


Proper  deployment  and  accurate  targeting  of  precision  guided  munitions  depend 
partly  on  knowledge  of  wind  velocity,  ceiling,  and  electro-optical  extinction 
in  the  target  area.  Knowledge  of  wind  in  the  target  area  is  crucial  for  the 
correct  placement  of  smoke  munitions,  since  wind  velocity  near  the  location 
where  the  artillery  is  deployed  may  differ  considerably  from  that  near  the 
target.  The  ability  to  accurately  estimate  atmospheric  characteristics  in  the 
target  area  also  should  permit  the  more  efficient  use  of  both  guided  and 
unguided  munitions  and  thereby  reduce  the  number  of  rounds  needed  to 

accomplish  a  specific  task. 

Wind  velocity,  ceiling,  and  volume  extinction  coefficient  may  be  computed  by 
methods  (developed  in  this  report)  that  use  only  data  currently  acquired  by  a 
remotely  piloted  vehicle  (RPV).  No  new  instrumentation  is  required;  the  input 
consists  of  flight  data  and  data  from  an  on-board  imaging  system.  Since  some 
of  the  principal  uses  of  RPV  now  and  in  the  near  future1  are  surveillance  of 
enemy-held  territory  and  target  detection  and  designation,  techniques 

developed  herein  may  be  used  to  describe  the  above  atmospheric  variables  in 

the  data  silent  region  near  the  target  while  tne  aircraft  performs  other 

missions  such  as  surveillance.  (See  Robinson*  for  a  general  description  of  a 
number  of  RPV  and  their  instrumentation  and  El  son*  for  information  on  the  RPV 
system  being  developed  for  the  Army.) 

The  data  for  the  first  algorithm  for  wind  velocity,  which  uses  along-wind 

information,  consists  of  heading,  airspeed,  and  ground  speed  (or  distance 

flown  and  time  to  fly  that  distance)  for  two  perpendicular  courses.  The 
second  wind  velocity  algorithm,  which  uses  crosswind  information,  requires 
heading,  airspeed,  and  drift  or  correction  angle  for  two  perpendicular 

courses.  Input  for  the  algorithm  for  computation  of  volume  extinction 
coefficient  consists  of  horizontal  flight  or  ground  distance,  altitude,  angles 
between  the  vertical  and  the  line  of  sight  (LOS)  to  the  radiating  surface,  and 
radiances  or  equivalent  voltages  from  the  radiating  surface  to  the  RPV  over 
two  separate  paths.  The  ceiling  algorithm  includes  horizontal  flight  or 

ground  distance,  altitude,  angle  between  the  flight  path  and  the  LOS  to  the 
cloudbase,  and  the  angle  between  the  flight  path  and  the  LOS  to  a  landmark 
vertically  below  the  viewed  cloudbase.  These  methods  are  embodied  in  computer 
programs  that  can  be  run  on  a  desktop  computer.  An  alternate  shorter  program 
is  presented  for  calculation  of  wind  velocity  when  both  headings  are  known. 
The  computer  codes,  in  BASIC,  are  shown  in  the  appendix. 


*M.  H.  Crowell,  1980,  A  Survey  of  Simulation  and  Test  Results  for  Assessing 
RPV  Performance  in  a  ~WBlC  Environment,  Pinal  Report  SPC  615,  prepared  by 
Systems  Planning  Corporation  for  PMi),  Tactical  Airborne  Remotely  Piloted 
Vehicle/Drone  Systems,  US  Army  Aviation  R4D  Command,  Contract  DAAK50-80- 
C-0011,  Saint  Louis,  MO,  50  pp 

*A.  Robinson,  1980,  "Battlefield  Reconnaissance:  Penetrating  the  Fog  of  War," 
Mil ita ry  Tech  and  Econ,  4 ( 1 5 ) : 3 3-4 2 

’B.  M.  El  son,  1980,  "Mini -RPV  Being  Developed  for  Army,"  Aviation  Week  and 
Space  Technology,  7  January  1980,  pp  2-7 
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2.  ALGORITHMS 

The  algorithms  described  in  this  report  are  mathematically  simple  and  easy  to 
understand.  The  manual  versions  of  these  procedures  require  only  a  simple 
calculator  or  trigonometric  tables,  graph  paper,  a  fine-scaled  ruler,  and  a 
pencil  and  paper.  A  desktop  computer  able  to  use  BASIC  is  sufficient  to 
perform  the  automated  versions.  The  operator  need  only  type  In  the  requested 
quantities. 

2.1  Wind  Velocity 

2.1.1  First  Technique:  Distance  and  Time  Input 

The  input  data  are  (1)  airspeed  in  meters  per  second  along  the  Initial  course 
and  the  second  perpendicular  path  (X  and  Y),  (2)  distance  in  meters  along  each 
path  and  the  associated  time  in  seconds  (Dx,  Tx  and  Dy,  Ty),  and  (3)  heading 

(direction  to)  of  each  course  (Dr,  and  Dr2).  To  obtain  ground  speed  (Xg  and 

Yg),  simply  divide  the  distances  by  the  appropriate  times  (Xg  =  Dx/Tx,  Yg  = 

Dy/Ty).  Subtracting  X  from  Xg  gives  the  difference  C;  similarly  Yg  -  Y  =  D. 

The  windspeed  (V)  is  computed  from  the  formula  for  the  hypotenuse  of  a  right 

triangle.  V  =  (C  +  D  )  ^  . 

The  computation  of  wind  direction  (Dir)  in  the  desktop  computer  version 
requires  values  of  C,  D,  Or!,  and  0r2.  If  either  Dri  or  0r2  is  missing  (input 
a  999  for  the  missing  value)  two  values  of  Dir  are  computed,  one  of  which  is 
correct.  The  correct  value  may  be  determined  with  the  aid  of  other 
information  (for  example,  a  synoptic  chart  can  indicate  which  of  the  two 

values  is  most  likely).  When  both  Dri  and  0r2  are  not  available,  a  message  is 
printed  saying  that  no  directions  were  given  or  computed.  A  further  condition 
for  the  computation  of  direction  is  whether  the  orientation  of  the  flight 
paths  is  "right"  or  "left."  In  the  context  of  this  report,  the  orientation  is 
determined  by  whether  the  Y  vector  is  to  the  right  or  left  of  the  X  vector 
when  facing  the  direction  of  flight  along  the  X  vector  (that  is,  toward 

Dri).  Numerically,  "right"  occurs  when  Dr2  >  Dri  (360°  added  to  Dr2  if  Dr2  < 
90°  and  270°  <  Dri  <  360°)  and  "left"  occurs  when  Drt  >  Dr2  (360°  added  to  Dr, 
if  Dr,  <  90°  and  270°  £  Dr2  <  360°). 

A  "flowchart"  (partly  in  plain  English)  can  provide  a  better  understanding  of 
the  intricacies  of  the  first  method  than  a  written  explanation  which  could  be 
tedious  and  somewhat  confusing  for  the  reader.  Such  a  chart  is  presented  in 
figure  1.  Figure  2  Illustrates  the  computation  of  wind  velocity  for  a  left 
orientation  when  C  >  0  and  D  <  0,  in  the  case  of  both  Dr,  and  Dr2  known. 

These  two  figures  should  be  used  together  to  gain  a  basic  understanding  of  the 

first  technique. 

The  present  form  of  this  algorithm  uses  distance  flown  and  time  to  fly  that 
distance  to  compute  ground  speed.  Simple  modifications  to  the  program  will 
permit  ground  speed  to  be  input  directly.  The  input  routine  would  need  a 
slight  modification,  and  the  simple  computation  of  ground  speed  would  be 
eliminated. 
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Neither  missing 


Dr2  missing 


Or2  *  Orj  -  90° 


Is  Dr?  <  0? 


Yes 


Dr2  «  Dr2  +  360° 


Is  0r2  >  360°? 


Yes 


Dr2  *  0r2  -  360 


Input  for  two  perpendicular  courses: 

(1)  Airspeeds  in  meters/second  (X,Y) 

(2)  Distances  in  meters  (Dx,  Dy) 
and  times  In  seconds  (T„,  Ty) 

(3)  Headings  (Drj,  Dr2) 

:l: z  =z 

Compute  ground  speed 
xg  1  DX/TX 
Yg  *  Dy/Ty 

~  _  r  ~ ~zz 

Compute  parallel  components 
of  the  wind 
C  =  X.  -  X 
0  =  Yg  -  Y 

i  . 

Compute  the  windspeed 

_ V  =  (C2  +  o2)1'2 _ 


Compute  angles  (P,  Q)  relating 
wind  and  flight  directions 
P  *  arctan  |C/D| 

Q  =  arctan  |D/C| 


Set  flag  for  C  and  D  with  the 
same  sign  (L  *  1)  or  opposite 
signs  (l  *  0) 


Is  Drj,  and/or  Dr2  unknown? 

] - 1 - 


Drj  missing 


Set  0 r2  =  Drj  +  90 

| 

Set  Drj  ■=  Dr2  +  90° 

zzj  . 

- ^ - 

Is  this  step  in  the 
first  or  second  pass? 

h“ 

Is  this  step  in  the 
first  or  second  pass? 

S5| 


Drj  »  Dr2  -  90° 


Is  Drj  <  0? 


es 


Is  Drj  >  360° 


es 


Drj  *  Drj  -  360 


Both  missing 


Print  “no  direction  given 
or  computed" 

Print  value  for  windspeed 


, .  Flowchart  of  the  first  technique  for  computation  of  wind  velocity. 
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0°  or  360° 


Figure  2.  Illustration  of  the  computation  of  wind  velocity  by  the  first 
algorithm  where  the  orientation  is  "left,"  C  >  0  and  D  <  G,  and 
both  Dr^  and  Dr2  are  known.  X  and  Y  are  airspeeds  and  Xg  and  Yg 
are  ground  speeds.  C  =  Xg  -  X  and  D  -  Yg  -  V.  V  is  the  wind  vector 
and  P  and  Q  are  computed  by  the  arctangents  of  the  absolute  values 
of  C/D  and  0/C,  respectively.  Dr^  and  Oro  are  the  directions 
toward  which  the  RPV  flies  along  the  X  and  Y  flight  paths, 
respectively.  In  this  example,  the  wind  is  blowing  from  a 

direction  slightly  less  than  90°. 
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2.1.2  Second  Technique:  Drift  or  Correction  Angle  Input 

The  input  data  are  (1)  airspeed  in  meters  per  second  along  the  two 
perpendicular  courses  (X  and  Y),  (2)  a  flag  to  indicate  whether  drift  (1=1) 
or  correction  (I  =  -1)  angles  are  utilized,  (3)  the  drift  or  correction  angles 
(A  and  8)  in  degrees,  and  (4)  the  heading  of  each  course  (Dri  and  0r2).  An 
angle  is  considered  positive  if  it  describes  an  arc  running  to  the  right  of 
the  flight  path  when  facing  in  the  direction  of  flight.  The  tangent  of  the 
angle  times  the  relevant  airspeed  gives  the  crosswind  component  for  each 
course.  C  =  X  tan  A,  D  =  Y  tan  B.  The  windspeed  (V)  is  computed  from  the 

j  2  i /2 

same  formula  as  in  the  first  algorithm;  that  is  V  =  (C  +  D  )  ' 

The  computation  of  wind  direction  (Dir)  requires  C,  D,  Dri,  and  Dr2  as 

input.  If  either  0rt  or  0r2  is  missing  (input  a  999  for  the  missing  value) 
two  values  of  Dir  are  computed,  one  of  which  is  correct.  When  both  Dri  and 
Dr2  are  missing,  a  message  is  printed  saying  that  no  directions  were  given  or 
computed.  The  orientation,  left  or  right,  is  determined  as  in  the  first 
method. 

A  "flowchart"  similar  to  that  of  figure  1  is  presented  in  figure  3,  but  for 

the  second  method.  Figure  4  illustrates  the  computation  of  wind  velocity  for 

a  right  orientation  where  C  >  0  and  D  <  0,  when  both  Dri  and  Dr2  are  known. 
These  two  figures  should  be  used  together  to  gain  a  basic  understanding  of  the 
second  technique. 

2.1.3  A  Shorter  Version 

A  shorter  version  of  the  computer  program  was  developed  that  has  about  two- 
thirds  the  number  of  statements  and  storage  requirement  as  the  program 
described  in  sections  2.1.1  and  2.1.2  of  this  report.  To  reduce  the  number  of 
statements,  it  was  assumed  that  Drt  and  Dr2  would  always  be  known.  If  either 
direction  is  unknown,  this  program  can  be  run  twice  with  the  unknown  direction 
=  the  known  direction  +90°.  (One  of  the  two  velocities  will  be  correct.) 
Consequently ,  all  statements  associated  with  the  extra  computation  required  to 
handle  the  cases  of  either  Dri  or  0r2  unknown  were  removed,  along  with  those 
activated  when  both  directions  were  missing.  A  listing  of  the  shorter  version 
is  in  the  appendix  along  with  the  listing  of  the  complete  version.  Figures  2 
and  4  illustrate  the  output  from  the  shorter  program. 

2.2  Volume  Extinction  Coefficient 

The  computer  programs  for  the  calculation  of  volume  extinction  coefficient 
require  the  input  of  either  (1)  horizontal  flight  or  ground  distance,  (2) 
altitude,  (3)  or  both  distance  and  altitude.  If  either  distance  or  altitude 
is  unknown,  then  (4)  the  angle  between  the  vertical  and  the  slant  path 
(technique  A)  or  (5)  the  angles  between  the  vertical  and  the  two  slant  paths 
(technique  8)  are  input.  Finally  (6)  radiances  or  (7)  equivalent  voltages  are 
entered  for  the  respective  views  of  the  radiating  surface(s). 

The  equations  are  derived  by  first  assuming  that  the  vertical  distribution  and 
amount  of  scatterers  and  absorbers  in  any  vertical  column  of  the  same  height 
are  constant  over  the  area  of  interest.  This  assumption  is  reasonable  to  a 
fair  accuracy  over  small  areas  of  the  order  of  a  few  tens  of  square  kilometers 
or  less,  not  in  the  Immediate  vicinity  of  atmospheric  "discontinuities"  such 
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START 


T 


Input  for  two  perpendicular  courses: 

(1)  Airspeed  In  ms"1  (X,  V) 

(2)  Whether  using  drift  (I  *  1)  or  correction  (I  *  -1)  anqle 

(3)  Angles  In  degrees  (A,  B) 

+  If  to  right  of  direction  of  flight 
-  If  to  left  of  direction  of  flight 

(4)  Headings  (Orj,  0r2) 


I 


A  =  I  x  A 
B  =  I  x  B 


Compute  perpendicular  components  of  the  wind 


C  *  X  tan  A 
D  *  Y  tan  B 


Compute  wlndspeed 
V  =  (C2  +  D2)1'2 


I 


Compute  angles  relating  wind  and 
flight  directions 
P  *  arctan  IC/Di 
Q  =  arctan  10/Cl 


Neither  missing 


0r2  missing 
[Set  0r2  -^Dr^  +  90° | 


Is  this  step  In  the 
first  or  second  pass? 


Second  ^ 


Set  flag  for  C  and  0  with  the  same  sign  (L  *  1) 
or  opposite  signs  (L  =  0) 


Dr2  =  Dr,  -  90' 


1 —  , 

I  "7<0 

Yes  j 


0r2  =  0r2  +  360” 


T 


Yes 


Is  Pr9  >  360* 

■  ir 


[Pr2  1  Pr2  -  360 


5- 


First 


No 


No 


1 

^Is  Drt  and/or  Pr2  unknown?  [_ 


Both  missing 


Orj  missing 
Set  Orj  =  0r2  +  90'  | 


Print  "No  direction  given 
or  computed" 

Print  value  for  windspeed 


I 


Is  this  step  In  the  1 _ _  First 

first  or  second  pass? 


Second 


Orj  =  0r2  -  90° 


I 


Is  Orj  <  0 


Tes  f 

[ Drt  »  Pr2  7  360'  [ 


I 


[  Is  Drx  >  360' 
Yes  X 


Orj  =  Orj  -  360' 


No 


No 


Figure  3.  Flowchart  of  the  second  technique  for  the  computation  of  wind  velocity. 
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Determine  whether  headwind 
or  tailwind 

(1)  If  Dri  <  90°  and  270°  <  Dr 2  <  360* 
then  Drj  =  Drj  +  360° 

If  Dr2  <  90°  and  270°  <  Drj  <  3bO* 
then  9rj  =  3rj  +  360° 

(2)  Tailwind  if: 

(a)  D  >  0  and  Drj  >  Dr2 

(b)  0  <  0  and  Dr^  <  Dr2 


Determine  whether  headwind 
or  tailwind 

(1) 

Same  as  for(T) 

(3) 

Tailwind  if: 

(a)  C  <  0  and  Drj 

>  Ilr2 

( b )  C  >  0  and  Drj 

<  Dr2 

Determine  orientation 
"Right"  if  Orj  <  Dr2 
"Left"  i f  Drj  >  Dr2 


Determine  orientation 
as  for  (a) 


Is 

orientation 

“right" 

and 

L  -  1 

(C  and  D  have  same 

sign)? 

Or 

Is 

orientation 

"left" 

and 

L  =  0  (C 

and  D  have 

O 

O 

c  signs)? 

Is  orientation  "right"  and 
L  =  1? 
or 

Is  orientation  “left"  and 
l  =  0? 


Dir  =  Dr,  ♦  PI 


I  Dir  -  Dr? 


[ Is  there  a  tailwind?! 


;  1 
|Dir  =  Di 

r  HW°  I 

I  Is  Dir 

>  3o0  ■?  t: 

|Dir  =  Dir  -  360“| 

1  -  --  k-J 

1  Is  Dir 

■  <  o? 

«  1 

ID i r  =  Di  - 

r  ♦  3b0°  | 

Is  the  input  value  o 

and  does  this  step 

f  Drj  or  Dr2  unknown 
end  the  first  pass? 

fes  , 

|  Return  to 

(*)  uhovel 

Print  values  of  windspood  and  direction 
(2  values  f  each  if  Drj  or  Dr2  unknown) 

i igure  3  (cont ) 


Figure  4.  Illustration  of  the  computation  of  wind  velocity  by  the  second 
algorithm  where  the  orientation  is  "right,"  C  >  0  and  D  <  0,  and 
both  Dr^  and  l)r2  are  known.  The  variables  shown  are  the  same  as  in 
figure  2  except  that  A  and  D  are  the  drift  angles  (-  -correction 
angles)  for  the  X  and  Y  flight  paths,  respectively.  The  crosswind 
components  are  C  =  X  tan  A  and  B  =  Y  tan  B.  In  this  example,  the 
wind  is  from  slightly  less  than  360°. 
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as  a  sharp  front.  Therefore,  between  a  given  flight  level  and  the  ground,  a 
change  in  total  mass  of  absorber  or  scatterer  roughly  is  deoendent  only  on  the 
difference  in  path  length. 

Taking  the  ratio  of  two  radiances  from  the  same  source  but  over  different  path 
lengths  and  using  Beer's  law,  we  have: 


-k'mZ*  -k'mZt 
R,/R,  =  R0e  /R0e 


where  R  =  radiance,  k1  =  mass  extinction  coefficient,  m  =  mass  of  absorber  and 
scatterer  per  unit  volume,  Z,  and  h  are  path  lengths,  and  the  subscripts  0, 
1,  2  refer  to  values  at  the  source  and  at  the  sensor  for  the  two  paths, 
respectively.  If  we  let  the  volume  extinction  coefficient  (k)  =  k'm,  and  we 
factor  out  R«,  we  have 


-kZj  -kZ1 
R ,/R ,  =  e  /e 

-kZ,  +  kZ, 
=  e 

k(Z,  -  Z,) 
=  e 


Taking  the  logarithm,  we  have: 

In  (R,/R, )  =  k  (Z,  -  Z,) 


and  for  k: 


k  =  In  (Rt/Rj )/(Z ,  -  Zt) 


If  the  RPV  can  fly  directly  over  the  radiating  surface,  we  can  use  technique  A 
(see  figures  5a  and  5b).  If  overflight  is  not  possible,  then  use  technique  B 
(figure  5c).  These  techniques  are  described  below;  the  reader  should  refer  to 
the  appropriate  figure  (5a,  b,  or  c)  in  the  following  descriptions. 

2.2.1  Technique  A:  Overflight  of  Target 

Here  we  assume  overflight  is  possible.  The  RPV  observes  two  surfaces  of 
closely  similar  properties  such  as  different  regions  of  a  lake  with  a  very 
nearly  uniform  surface  temperature,  or  It  looks  at  the  same  surface,  once 
vertically  and  again  along  a  slant  path. 
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X 


Figure  5a.  The  geometry  for  finding  the  volume  extinction  coefficient 

(k  =  k'm)  by  viewing  the  same  target,  (k1  =  mass  extinction 
coefficient,  m  =  mass  of  absorber  and  scatter  per  unit  distance.) 
x  =  ground  or  horizontal  distance,  8  =  angle  between  vertical  and 
line  of  sight  to  target,  and  Z  and  Z2  are  the  vertical  and  slant 
paths. 


Figure  5b.  Finding  k  by  viewing  two  closely  similar  targets.  Variables  are 
as  in  5a. 
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Figure  5c.  Finding  k  by  viewing  the  same  target  along  different  slant  paths 
(different  anales,  both  >  0°).  Here  Z  =  altitude,  Zj  and  Z2  ~ 
slant  paths,  R^  and  R2  =  respecti ve  radiances ,  8j  and  e2  =  angles 
between  the  line  of  sight  to  the  target  and  the  vertical,  and  Xj 
X2  are  horizontal  distances  along  the  ground  or  flight  path. 
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Referring  to  figures  5a  and  b,  we  see  that  X  =  distance  along  the  ground,  Z  = 
Zi  =  altitude  (vertical  path).  It  =  slant  path,  and  s  =  angle  between  the 
vertical  and  the  slant  path. 

a.  If  X  and  Z  are  known  but  not  s,  then 

B  =  arctangent  (X/Z) 
and  It  ~  Z/cos  B. 

b.  If  Z  and  b  are  known  but  not  X,  then  as  in  a 

h  -  Z/cos  8. 

c.  If  X  and  8  are  known  but  not  Z,  then 

Z2  =  X/sin  8 
and  Z  =  X/tan  6 

Using  equation  (1)  and  the  above  geometrical  relationships,  we  have  for  the 
first  two  cases  (a,  b) : 


k  =  In  (R./RJ/ZU  -  1/cos  s) 
and  for  the  third  case  (c) 


(Z) 


k  =  In  (R,/R. )/X( 1/tan  8  -  1/sin  b) 


(3) 


2.2.2  Technique  B:  Standoff  from  Target 

Here  we  assume  overflight  is  not  possible.  The  RPV  views  the  same  surface 
from  two  different  angles  (neither  path  vertical  as  in  figure  5c)  or  views  two 
closely  similar  surfaces  (no  figure  shown). 

Referring  to  figure  5c,  we  see  that  X,  and  X,  =  horizontal  distances  along  the 
flight  path  or  the  ground,  Z  =  altitude,  Zi  and  h  are  slant  paths,  and  Bi  and 
8i  are  the  respective  angles  between  h  and  Z»  and  the  vertical.  Therefore: 

d.  If  Xi,  Xi  and  Z  are  known  but  not  8«,  B»,  then 

8i  =  arctangent  (Xi/Z),  6j  =  arctangent  (Xt/Z) 
and  Zi  =  Z/cos  8i,  h  -  Z/cos  8» 

e.  If  Z  and  Bi,  62  are  known  but  not  Xi ,  X2,  then  as  in  d 

Z»  =  Z/cos  81 ,  Z  =  Z/cos  6 2 
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f.  If  X2,  X2  and  Bi,  6 2  are  known*  but  not  Z,  then 
Zi  =  Xi/sfn  Si,  Z2  -  X2/sin  b2 
and  Z  =  Xi/tan  Bi  =  X2/tan  62 

Using  equation  (1)  and  the  above  geometrical  relationships,  we  have  for  cases 
d  and  e: 


k  =  In  (Rj/R 1 )/Z( 1/cos  0 , -  1/cos  B2) 


(4) 


and  for  case  f: 


k  =  in  (R2/Rj )/(X, /sin  B)  -  X2/sin  b2) 


(5) 


However,  for  fit  -►  0  (see  figure  5c)  equation  (5)  fails.  To  avoid  this 
problem,  replace  X,/sin  81  with  Z  (=  X2/tan  Bi)  for  small  values  of  81  (  = 
several  degrees  «*$s).  Therefore,  in  place  of  equation  (5)  use: 


k  =  In  (R2/Rl )/Xj (1/tan  b2  -  1/sin  b2) 


(6) 


Mote  that  equations  (6)  and  (3)  are  the  same,  where  X2  =  X  and  s 2  =  0  {compare 
figures  5a  and  5c). 


2.2.3  Computation  of  Radiance  from  Voltage 

Modifications  of  the  two  techniques  were  developed  where  radiance  is  not  input 
directly,  but  is  computed  from  voltages.  Commonly,  sensors  use  devices  that 
transform  received  energy  into  voltages,  which  in  turn  are  converted  into 

radiances  by  the  means  of  some  algorithm.  Although  the  form  of  the  conversion 
algorithm  may  differ  from  one  sensor  to  another  (for  example,  linear  or 

quadratic),  a  simple  linear  form  was  used  here  only  to  indicate  how  such 
equations  would  fit  into  the  techniques  described  in  this  report.  The 

equations  have  the  form  R  =  a  +  bV  where  R  =  radiance,  V  =  voltage,  and  a  and 

b  are  constants  determined  empirically  during  calibration.  It  is  assumed  that 
one  voltage  is  produced  for  each  of  the  two  views  of  the  radiating  surface  or 
target,  yielding  the  two  required  radiances. 


*These  computations  may  be  performed  if  either  b«  or  Bi  is  known.  For 
example.  If  Xi  and  Bi  are  known,  we  can  compute  Zi  and  Z.  Having  7.  and  X2,  we 
can  then  calculate  b2. 


2.2.4 
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Computer  Programs 

Four  "flowcharts"  were  constructed  to  enable  the  reader  to  understand  the 
computer  programs  more  easily  and  to  avoid  the  possible  confusion  an  entirely 
written  explanation  would  cause.  The  reader  should  refer  to  figure  5  when 
viewing  these  flowcharts.  Figure  6  shows  the  chart  for  technique  A,  where  the 
RPV  can  overfly  the  target  (=  radiating  surface)  and  radiances  are  input. 

Figure  7  presents  the  chart  for  technique  A  when  voltages  are  input.  Figure  7 

differs  from  figure  6  in  the  substitution  of  the  voltage  input  and  conversion 
algorithm  in  place  of  the  radiance  input  statement.  Figure  8  has  the 
flowchart  for  technique  8,  where  the  RPV  cannot  overfly  the  target  and 
radiances  are  input.  Figure  9  illustrates  the  difference  in  the  flowchart  for 
technique  B  when  voltage  input  replaces  radiance  input;  the  entire  flowchart 
is  not  shown.  The  programs  for  both  techniques  are  presented  in  the  appendix. 

2.3  Ceiling 

Ceiling  (c)  may  be  computed  by  using  simple  geometry  and  data  from  an  RPV 
carrying  a  movable  sensor  active  in  any  imaging  wavelength  region.  The 
required  input  includes  (1)  upward  elevation  angle  which  is  the  angle  between 
the  flight  path  and  the  LOS  to  cloudbase,  and  either  (2)  horizontal  flight  or 
ground  distance,  (3)  altitude,  or  (4)  both  altitude  and  distance.  If  only 

altitude  or  distance  is  known,  then  input  (5)  the  depression  angle  which  is 

the  angle  between  the  flight  path  and  the  LOS  to  a  landmark  vertically  below 
the  view  of  the  cloudbase. 

Referring  to  figure  10,  we  have  for  the  ceiling 

c  =  h  +  Z 

where  Z  =  altitude  and  h  =  vertical  distance  from  flight  level  to  cloudbase. 
Furthermore,  h  =  Xtan  b,  where  X  =  horizontal  flight  or  ground  distance  and  b 
=  upward  elevation  angle.  Substituting  for  h,  we  have  for  the  case  of  X  and  Z 
known: 


c  =  Z  +  Xtan  6  .  (7) 

If  X  is  known  but  not  Z,  we  have  for  Z 

Z  =  Xtan  a  , 

where  a  =  depression  angle.  Using  equation  (7)  and  substituting  for  7,  we 
have  for  c 


c  =  Xtan  a  +  Xtan  8 
=  X  (tan  a  +  tan  e)  . 
If  Z  is  known  but  not  X,  we  have  for  X 

X  =  Z  ctn  a 
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Figure  7.  Flowchart  of  technique  A  for  estimation  of  volume  evtlnctljn  coefficient  using  voltage  Input 


Figure  ft.  Flowchart  of  technique  P  for  estimation  of  volume  extinction  coefficient. 
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Using  equation  (7)  and  substituting  for  X,  we  have  for  c 


c  =  Z  +  Z  ctn  a  tan  8 
=  Z  (1  +  tan  s/tan  a)  . 

It  has  been  assumed  that  the  cloudbase  can  be  observed  directly  above  the 
landmark  from  the  RPV  (see  figure  10).  The  computation  of  c  becomes  less 
accurate  as  the  cloudbase  to  landmark  path  departs  from  vertical,  although  a 
departure  of  only  a  few  degrees  is  not  significant.  Furthermore,  it  is 
assumed  that  the  distance  X  is  the  same  for  both  views,  to  the  cloudbase  and 
to  the  landmark.  The  value  of  c  will  depart  from  the  real  value  as  the 
difference  in  X  for  the  two  views.  This  problem  could  be  solved  if  a  side 
looking  sensor  was  used  or  the  RPV  flew  at  a  low  speed  and  X  was  large.  For 
example,  relatively  little  degradation  will  occur  if  the  speed  of  the  RPV  =  20 
ms-1,  X  =  4000  m,  and  the  sensor  viewed  both  scenes  within  2  s. 

The  program  for  calculation  of  ceiling  is  presented  as  a  "flowchart"  in  figure 
11.  The  reader  also  should  refer  to  figure  10  as  an  aid  to  understanding  the 
flowchart.  The  computer  code  for  this  program  is  presented  in  the  appendix. 

3.  SAMPLE  COMPUTER  RUNS 

3.1  Wind  Velocity 

Four  examples,  two  for  each  algorithm,  are  presented  in  this  section  to  better 
demonstrate  the  computation  of  wind  velocity  by  the  two  methods.  For  each 
example,  a  table  shows  the  calculations  required  and  an  accompanying  graph 
shows  the  graphical  solution.  Each  set  of  one  table  and  one  graph  is 
presented  in  the  form  of  one  figure  for  ease  of  understanding  (figures  12 
through  15).  Although  more  than  four  situations  exist  (for  example,  "right” 
orientation  for  both  C  and  D  >  0  and  both  Dri  and  Dr2  known),  to  include  them 
all  for  both  methods  would  make  this  report  unnecessarily  large  and  tedious. 
Table  1  gives  the  computer  output  for  these  four  examples. 

Figures  12  and  13  present  solutions  for  the  first  algorilnm.  In  ffgive  12, 
C  >  0  and  D  <  0,  and  the  orientation  is  not  known  sines  only  Dr2  is  given. 
Two  solutions  are  computed,  one  of  which  is  correct.  Figure  13  has  a  "left" 
orientation  where  C  <  0  and  D  >  0.  Figure  2  and  section  2.1.1  describe  the 
relevant  variables. 

Figures  14  and  15  show  solutions  for  the  second  algorithm.  Figure  14  has  an 
orientation  that  is  "left"  and  C  and  D  are  negative.  In  figure  15,  the 
orientation  is  unknown  since  only  On  is  known,  and  C  and  D  are  both 
positive.  Figure  4  and  section  2.1.2  describe  the  relevant  variables. 

3.2  Extinction  Coefficient  and  Ceiling 

A  series  of  runs  of  the  ceiling  and  the  extinction  coefficient  (radiance 
input)  programs  were  made  to  illustrate  the  algorithms.  Sample  output  from 
the  extinction  programs  that  input  voltage  are  not  shown  here  because  they 
essentially  repeat  the  results  of  the  radiance  versions  and  a  lengthy  series 
of  examples  would  be  unnecessarily  tedious.  Tables  2  and  3  show  computations 
for  both  techniques  in  which  both  manual  and  computer-generated  values  are 
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shown.  Input  values  are  listed  with  the  appropriate  equations.  The  radiance 
values  have  no  specific  units  since  they  could  have  any  standard  (or 
nonstandard)  units  without  affecting  the  results.  In  any  case,  the  ratio  of 
the  radiances  is  dimensionless.  Also,  the  values  of  k  in  these  tables  were 
computed  for  comparison  and  technique  demonstration  purposes,  assuming  perfect 
input.  In  the  real  world,  the  last  two  or  three  digits  to  the  right  of  the 
decimal  point  probably  would  be  meaningless. 

Ceiling  computations  are  shown  in  table  4  in  a  format  similar  to  that  of 
tables' 2  and  3.  Input  values  are  listed  along  with  the  appropriate  equations 
and  both  manual  and  computer-generated  values  are  shown.  Finally,  table  5 
presents  samples  of  output  from  the  computer  programs  used  to  generate  the 
values  in  tables  2  through  4,  and  output  from  the  "voltage"  versions  of  the 
extinction  programs. 


Figure  10.  Geometry  for  estimation  of  ceiling  (  c  =  h  +  z)  where  z  = 

altitude,  h  =  height  of  cloudbase  above  flight  level,  x  =  ground 
distance,  and  a  and  e  are  angles  relating  z  and  h  to  x, 
respectively. 
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Figure  11.  Flowchart  for  the  estimation  of  celling. 
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<  -  o9  ..i, '  .  j 

r  *  /6  Mb'*  •{  -  iij- 


A  =  I  x  A  =  -3° 

B  -  I  x  B  =  -10° 

C  =  X  tan  A  -  -3.1 

D  -  Y  tan  8  =  -13 

V  =  (C2  +  02)1/2  =  14.3  ms' 


P  =  arctan  |  C/D| 
0  -  arctan  |D/C| 

Dir  =  0r:  +  P  =  104.6° 


0  or  360° 


Pr£ 


Figure  14.  Use  of  method  2  to  compute  wind 
known.  Variables  are  explained 
3oth  C  and  D  are  negative.  T1 
"left." 
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Y  =  50  ms 


B  =  -25 


Dr2  »  999°  (unknown) 


A  =  I  x  A  =  15° 

8  =  I  x  B  =  25“ 

C  =  X  tan  A  1  13.4  ms"1 

0  =  Y  tan  B  *  23.3  ms"1 

V  =  (C2  +  02)1/2  =  26.9  ms'1 

P  =  arctan  |C/0|  =  29.9° 

Q  =  arctan  | 0/C[  =  60.1° 

0r2  =  Dr1  t  90°  =  335°  or  155° 

(1)  Dir  =  Orj  -  P  =  215.1° 

(2)  Dir  =  0rL  +  P  ♦  180  =  454.9  -  360  =  94.9° 


180° 


Figure  15.  Use  of  method  2  to  compute  wind  velocity.  Variables  are  explained 
in  figure  4  and  second  technique.  Only  Orj  is  known,  and  C  and  D 
are  both  positive.  More  than  one  value  of  some  variables  is  shown 
because  Dr2  may  be  155°  or  335°.  Two  wind  velocities  (V)  are 
computed. 
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TABLE  1.  COMPUTER  OUTPUT  FOR  THE  FOUR  EXAMPLES  OF  FIGURES  12  THROUGH  15. 

AS  NOTED  IN  THOSE  FIGURES,  TWO  ANSWERS  ARE  GIVEN  WHEN  ONE  HEADING 
IS  UNKNOWN.  THE  MOST  PROBABLE  OF  THE  TWO  WIND  DIRECTIONS 
IS  DETERMINED  WITH  THE  AID  OF  OTHER  DATA  (FOR  EXAMPLE,  A 
SYNOPTIC  CHART). 


5  WIND  VELOCITY 


Windspeed  =  39.1  ms' 
Windspeed  =  39.1  ms' 


Wind  direction  =  245.2° 
Wind  direction  *  324.8° 


ONE  ANSWER  IS  CORRECT 


6  WIND  VELOCITY 


Windspeed  =  11.7  ms'* 


7  WIND  VELOCITY 


Windspeed  *  14.2  ms'* 


8  WIND  VELOCITY 


Windspeed  =  26.9  ms'1 
Windspeed  =  26.9  ms'* 


Wind  direction  =  1.6° 


Wind  direction  =  104.7° 


Wind  direction  *  215.1° 


Wind  direction  =  94.9 


ONE  ANSWER  IS  CORRECT 


TABLE  2.  SAMPLE  COMPUTATIONS  OF  VOLUME  EXTINCTION  COEFFICIENT  USING  TECHNIQUE 
A  WITH  RADIANCE  INPUT,  k  *  VOLUME  EXTINCTION  COEFFICIENT,  X  = 
HORIZONTAL  DISTANCE,  Z  =  ALTITUDE,  8  =  ANGLE  BETWEEN  VERTICAL  AND 
LOS  TO  TARGET,  AND  R  AND  Rj  =  RADIANCES  FROM  THE  TARGET  WITH 
VERTICAL  AND  SLANT  VIEWING,  RESPECTIVELY.  FOR  THESE  EXAMPLES 
R  >  100  AND  Rj  *  25  UNITS;  THEREFORE,  Rj/R  *  0.25  and  In  (Rj/R)  = 
-1.3863.  NOTE  THAT  HERE  AND  IN  THE  COMPUTER  PROGRAMS  Rj  *  R2  OF 
EQUATION  (3)  AND  R  =  Rj  OF  EQUATION  (3). 

Known  Variables  Equations  and  "Manual11  Values  Computer  Values 

a)  X  =  1000  m  *  1.0  km  e  -  arctan  (X/Z) 

Z  *  200  m  =  0.2  km  =  78.69° 

k  =  In  R|/R/Z(l  -  1/cos  8) 

=  -1.3863/0.2(1  -  5.0990) 

=  1.6910  km'*  1.6910  km'1 


b)  Z  =  300  m  =  0.3  km  k  =  -1.3863/Z(1  -  1/cos  8) 

8  =  65°  =■  -1.3863/0.3(1  -  2.3662) 

=  3.3824  km'1 

c)  X  *  600  m  =  0.6  km  k  =  -1 .3863/XI 1/tan  8  -  1/sin  e) 

8  =  50°  =  -1.3863/0.6(0.8391  -  1.3054) 

=  4.9550  km'* 


3.3824  km'* 


4.9549  km'* 
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TARLE  3.  SAMPLE  COMPUTATIONS  OF  VOLUME  EXTINCTION  COEFFICIENT  USING  TECHNIQUE  V. 

WITH  RADIANCE  INPUT.  k  =  VOLUME  EXTINCTION  COEFFICIENT,  Xj  AND  X2  = 
HORIZONTAL  DISTANCE,  7.  =  ALTITUDE,  Sj  AMD  B2  =  ANGLES  BETWEEN  VERTICAL 
AMO  LOS  TO  TARGET,  AND  R,  AND  R?  *  RADIANCE  FROM  THE  TARGET  ALONG  THE  TWO 
SLANT  PATHS.  FOR  THESE  EXAMPLES.  Rj  =  100  AND  R2  =  30  UNITS;  THEREFORE, 
R?/Rj  =  0.30  AND  In  (R,/R2)  =  -1.2040. 


Known  _Variabl_es_ 
a)  Xj  -  300  n  =  0.3  kn 
X,  =  1500  m  =  0.5  kn 
Z  =  500  r,i  =  0.5  km 


h)  Z  -  500  m  =  0.5  km 

ij  =  20° 

a2  =  70° 


Equations  and_^'ManuaV^^alues  Computer  Values 

8j  -  arctanqent  (Xj/7)  -  30.96° 

32  =  arctanqent  (X?/7)  =  71.57° 

k  =  In  (Rp/Rj )  'Z(  1  /cos  9j  -  1/cos  ,;) 

=  -1.2040/0.5(1.661  -  3.1631) 

-•  1.205H  km'1  1.2063  km'1 

k  =  -1.2040/0.5(1.0642  -  2.0235) 

=  1.2949  km'1  1.2944  cai'1 


c)  Xj  =  200  m  =  D.2  km 


X2 

!1 

■’ 2 


2000  m  =  2.0  km 
15° 

69.63° 


k  =  1.2040./;<j/sln  ij  -  Wsin  ii2) 

’  -1.2040/(0.7/27  -  2.1347) 

=  0.5540  km'1  0.5539  km'1 


3)  Xj  =  40  m  =  0.04  km 
X,  =  1500  m  =  1.5  km 
U  *  3° 
a,  -  63.03° 


k  =  -1  . 2040/ X 2 ( 1/tan  <2  -  1/sin  <2) 

-  -1.2040/1.5(0.5089  -  1.1220) 

=  1.3091  km'1  1.3090  km'1 


TABLE  4.  SAMPLE  COMPUTATION  OF  CEILING,  c  =  CEILING  (HEIGHT  OF  CLOUOUASE ) , 

z  =  altitude,  h  =  height  oe  clouobase  above  flight  path,  x  * 

HORIZONTAL  0! STANCE,  3  =  UPWARD  ELEVATION  ANGLE  BETWEEN  FLIGHT  RATH 
ANO  LOS  TO  CLOUDBASE,  ANO  i  =  DEPRESSION  ANGLE  BETWEEN  ri. IGHT  PATH 
AMO  LOS  TO  LANDMARK  (VERTICALLY  BElOW  VIEW  OE  ClOiJDBASE).  COMPUTE  I 
VALUES  OF  c  ARE  TO  THE  NEAREST  METER. 

Known  Variables  Equations  and  "Manual"  Values  Compute  Values 


a) 

X 

-  700  lA 

c  =  Z  ♦  Xtan  a 

7 

=  150  m 

=  15.)  *  233 

\ 

=  22’ 

=  433  n 

433  m 

t>) 

i 

-  300  m 

c  =  Z  ( 1  +  tan  a /tan  i ) 

3 

■  14° 

-  303(1  ♦  0. 24933/0. l’:>33) 

l 

=  10° 

=  724  m 

724  m 

c) 

X 

=  1000  m 

c  -  X( tan  i  *  tan  a ) 

=  20° 

=  1000(0.5/735  +  0.36397) 

=  30° 

-  941  m 

941  m 
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TABLE  5.  SAMPLES  OF  OUTPUT  FROM  THE  COMPUTER  PROGRAMS  FOR  CEILING  AND  VOLUME 


EXTINCTION  COEFFICIENT.  ONLY  THE  ANSWER  (FOR  EXAMPLE,  Ceiling  = 
100  m)  IS  PRINTED  BY  THE  COMPUTER.  SINCE  THE  DESKTOP  COMPUTER 
CANNOT  PRINT  SUPERSCRIPTS,  km'1  IS  PRINTED AS  1/km.  _ 


1.  Ceiling  given  altitude  (Z)  =  150  m,  horizontal  distance  (X)  =  700  m,  and 
elevation  angle  (?)  =  22°. 

Ceiling  =  433  m 

2.  Volume  extinction  coefficient  using  one  slant  path  given  Z  =  200  m,  X  = 

450  m,  and  radiances  (R  and  R^)  =  110  and  25  units. 

Volume  extinction  coefficient  =  5.0663  1/km 

3.  Volume  extinction  coefficient  using  two  slant  paths  given  horizontal 
distances  (Xj  and  X2 )  =  200  and  1800  m,  one  of  the  two  angles  between  the 
LOS  to  the  target  and  the  vertical  (6^)  =  18°,  and  radiances  (3^  and  = 
100  and  20  units. 

Volume  extinction  coefficient  =  1.2823  l/km 

4.  Volume  extinction  coefficient  using  one  slant  path,  and  using  voltages  as 
input  given  Z  =  200  m,  X  =  450  m,  voltage  coefficients  (A  and  B)  =  2.2  and 
11.0,  and  voltages  (Vj  and  V2)  =  8.0  and  2.0  v. 

Volume  extinction  coefficient  =  4.4989  l/km 

5.  Volume  extinction  coefficient  using  two  slant  paths,  and  using  voltages  as 
input  given  Z  =  200  m,  X^  and  X2  =  200  and  450  m,  8^  =  45°,  A  and  B  =  2.2 
and  11.0,  and  V^  and  =  8.0  and  2.0  v. 

Volume  extinction  coefficient  =  6.2771  l/km 


4.  CONCLUSION 

Useful  tools  for  silent  area  analysis  have  been  developed  in  the  form  of 
simple  methods  for  the  computation  of  wind  velocity,  ceiling,  and  volume 
extinction  coefficient.  These  algorithms  use  information  already  gathered  by 
an  RPV  of  the  type  being  developed  for  the  Army;  no  new  instrumentation  is 
required.  Windspeed  should  be  accurate  to  several  tenths  of  a  meter  per 
second  and  wind  direction  to  less  than  a  degree,  assuming  that  the  input  is 
"perfectly"  accurate.  Similarly,  ceiling  should  be  correct  to  about  1  or  2 
percent  and  extinction  coefficient  to  about  10  percent.  However,  under 
operational  conditions,  the  accuracy  of  the  input  data  probably  would 
determine  the  accuracy  of  the  output. 

Computations  may  be  performed  via  a  desktop  computer  able  to  use  the  BASIC 
computer  language,  or  by  the  use  of  a  hand-held  calculator,  a  fine-scaled 
ruler,  and  graph  paper.  The  former  technique  only  requires  the  operator  to 
input  numbers  that  are  specifically  requested;  the  latter  manual  technique 
requires  some  knowledge  of  the  situation. 
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APPENDIX 

COMPUTER  CODES  IN  BASIC  FOR  THE  RPV  PROGRAMS 

(A)  Is  the  full  version  and  (B)  Is  the  shortened  version  of  the  wind  velocity 
program;  (C)  Is  the  celling  program;  (D)  and  (E)  are  the  extinction  coeffi¬ 
cient  programs  with  one  and  two  slant  paths,  respectively;  (F)  and  (G)  are  the 
extinction  coefficient  programs  using  voltage  Inputs,  with  one  and  two  slant 
paths,  respectively. 


10  ! 

20  1  COMPUTE  UIND  VELOCITY  USING  RPV  PROGRAM 
30  ! 

40  A*=“Uind  Speed  *  “ 

50  B$*"Uind  Direction  *  “ 

60  D*«“  m/s “ 

70  E*“"  Degrees" 

80  F*=“UIND  VELOCITY" 

80  N-0 

100  H=  0 

110  INPUT  “Input  L“1  to  compute  V  using  Method  1  or  L*2  for  Method  2“ ,L 
120  IF  L=2  THEN  GOTO  Meth 

130  INPUT  “Input  airspeeds  in  m/s  <X  and  Y>“,X,Y 

140  INPUT  "Input  distances(m)  and  tieesfs)  < Dx ,  Tx . Dy , Ty > " , Dx , Tx , Dy , Ty 

150  INPUT  “Input  headings  in  degrees  <Dr1  and  Dr  2 , Dr  1 , Dr2 

160  PRINT  USING  Heading;F* 

170  Heading;  IMAGE  15X,13A,2/ 

180  GOTO  Anda 

190  Meth;  INPUT  “Input  airspeeds  in  m/s  (X  and  Y)“,X,Y 

200  INPUT  "Input  whether  drift  < I = 1 >  or  correction  <I*-t)  angle", I 

210  INPUT  “Input  angles  in  degrees  <A  and  B)‘,A,B 

220  INPUT  "Input  headings  in  degrees  <Dr1  and  Dr2 >“ , Dr  1 , Dr2 

230  A*I*A 

240  B-I*B 

250  PRINT  USING  HeadingjF* 

260  GOTO  Otra 

270  Anda;  CALL  Compl < X, Y, Ox , Tx , Oy , Ty , Dr  1 , Dr2, V, Di r , N > 

280  GOTO  Printer 

290  Otra;  CALL  Comp2< X , Y . A , B , Dr  1 , Dr 2 , U , V , Di r , I , M > 

300  Printer;  PRINT  USING  Ti t le ; A* , V, D* , B», Dir , E* 

310  Title;  IMAGE  5X, 1 3A, ODD . D, 4A , SX, t 7A, DDO . D, 8A, / 

320  IF  Nxt  THEN  GOTO  Anda 
33u  IF  M=1  THEN  GOTO  Otra 

340  IF  <N=2>  OR  <M-2>  THEN  PRINT  "  ONE  ANSWER  IS  CORRECT" 

350  PRINT  USING  Out 
360  Out ;  IMAGE  5X,6/ 

370  END 
380  ! 

390  ! 

400  1  SUBPROGRAMS 

41  0  1 

420  ! 

430  SUB  Compl <  X,  Y,  Dx ,  Tx ,  By ,  Ty  .  Dr  1  ,  Dr2,  V  .  Dir  ,  N  > 

440  I 

450  Dir*0 

460  DEG 

470  L«0 

480  Xg=Dx/Tx 

490  Yg«0y/Ty 

500  C-Xg-X 

510  D»Yg-Y 

520  V=<C-2+D*2)A  .5 

530  IF  C-0  THEN  C- . 001 

540  IF  D«0  THEN  D-.001 

550  ! 

560  !  Computation  of  Direction 
570  ' 

580  P^ATN<A8S<C/D>> 

590  0»ATN<  ABS<  D/C  >  > 

600  IT  <C>0)  AND  ( D  >  0  i  OR  <C<0>  AND  < D< 0  >  THEN  L*1 
6 1 0  Oron=Dr 1 
620  DrtO’Dr2 

630  IF  Drt*999  THEN  GOTO  Second  1  Is  Drl  missing1 

640  IF  Dr2»999  THEN  Drto=Dr1+90  !  Is  Dr2  missing1 

->50  IF  <  Dr  2*999  >  AND  >N»1)  THFH  Drto«Dr1-90 
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660  IF  DrtoiG  'l  HEN  OrtC'-iiilu  KihU 

670  IK  Dr  to >36  0  I'HFN  Dr  F  c.-pr  t  c- .'6  0 

680  IF  (t'i'on<  =  '*(i)  AND  <  Dr  t  c- >=2  / 0  >  AND  (.  Dr  to  •  =  36  0  : 

690  IF  <Drto<=90>  AND  <Dror.>*2?0>  AND  <Dron<=36o- 

'00  Dir=Dr1-G  ‘  H 

710  IF  <  Dr to<  Dr on  >  AND  CL=0>  THEN  Dir=Drt*Q 

720  IF  <  Drto>Dr  or.  >  AND  <1=1  >  THEN  Oir=Dr1+Q 

730  IF  C>0  THEN  Dir=Dir+1S0  ’  T 

740  GOTO  Direct 

^50  Second :  IF  Dr2=999  THEN  GOTO  AH  1  h 

760  Dron=Dr 2+90  j 

770  IF  N~ !  THEN  I>ron-Dr2-90 

780  IF  Dront  O  THEN  Dron-OrorN  >60 

.•90  IF  Dron>360  THEN  Dron=Dr on -36 0 

8u0  IF  <  Dron<  =9  0  >  AND  < Drto>=2?0  >  AND  < Or t o< = ’AO  • 

810  IF  <  Dr  to<  =9  0  )  AND  <  Dr  or.  >--270  >  AND  <  0rer.<»3eC  / 

020  Dir =Dr2+P  1  H 

830  IF  >;  Dr  t  o<  Dr  on  >  AND  <l.*0  •  THFn  Dir«0r7-P 

840  IF  <  Drto  >Drr.n  )  AND  C  t  —  1  >  THFN  D  1  r  =Dr  7-p 

850  IF  D  '  0  THEN  Dir=Dir  +  1:s0  1  T 

86  0  GOTO  D i r ect 


370  Alt:  -RINT  "NO  DIRECTION  GIVEN  OR  COMPUTED" 

380  Direct:  IF  0ir>360  THEN  Dir=0ir-36u 

390  IF  D i r  <  0  THEN  Dir=Dir  +  360 

900  IE  <  Dr  1 =999  >  OR  <  Dr  2  =  999  >  THEN  N  =  N+1 

910  SUBEND 

92  0  ! 

970  > 

94 o  SUB  Comp2< ,  Y ,  A ,  6 ,  Dr  1  . Dr  2, U . V  .Dir . 1  .Mi 

9  5  0  ! 

960  Dir=0 

97  0  K=  0 

98  0  U  =  0 
990  DEG 

1  00  0  C=X*TAN<  A  ) 

1010  D  =  Y*T  ANt  6  ) 

1020  V-f C*24D‘2>* .f 
1030  IF  C  =  0  THEN  C=.GG1 
1040  IE  D  =  G  THEN  D=.001 
105  0  ! 

1060  !  Computation  of  Direction 
107  0  ! 

030  P=ATN<  ABS<  C/D  )  ) 

1 090  Q=ATN<  A8S<  D/C  >  > 

1100  IF  <C>0>  AND  <D>0>  OF  >.C<0>  AND  <  D<  0  >  THEN  k  = 
‘110  Dron=Dr  t 
1120  Dr to=Dr2 

1130  IF  Dr  1-999  THEN  GOTO  Second  1 

1140  IF  Dr 2=999  THEN  Drto=0r1+90 

1150  IF  <Dr2=999)  AND  <M=t>  THEN  Drto-Drt-90 

1160  IE  Dr to<  0  THEN  Drto=0rto  +  360 

: 1 70  IE  Drto  >360  THEN  Drto=0rto-36 0 

1180  IF  <Dron<-90>  AND  <Drto>=270>  AND  <Drto<=360> 
1190  IF  < Drto<  =90  )  AND  < Dr on >=27 0 >  AND  <Dron<«360> 
1200  IF  <  D  >  0  )  AND  <0ron>0rto,’  THEN  U=t  ! 

1210  IF  <  D<  0  )  AND  <Dron<Drto)  THEN  U=1  1 

1220  Dir=Dr1-P  ! 

1230  IF  <  Dron>Drto  >  AND  <K  =  1>  THEN  Dir=Dr1+P 
1240  IF  < Dron<  Dr  to  7  AND  <K  =  0>  THEN  Dir  =  Drt+P 
1250  IF  U=1  THEN  Dir=*Dir  +  )80  ! 

1260  GOTO  Direct 

1270  Second:  IF  Dr2=999  THEN  GOTO  Alt  ! 

1280  Dron=Dr2+90  ! 

1290  IF  M=1  THEN  Dron*Dr2-90 
1300  IE  Dron<0  THEN  Dr on=Dr on+36 0 
'710  IE  Dron>360  THEN  Dr  on -Dr  or. -360 


THEN  L‘r  on=i>r  ■; r.+  in-C 
THEN  Drt<w=Dr  t  o+  76  0 
e  a  d ...  i  r,  d 


j  i  1 ...  i  i  i.J 

re  C'r  1  jj,  Dr 2  missing 
Dr  i  i-:  »i  i;ir,q 


IF'EN  Dt  or>=Dror  +  360 
THEN  Drto=0rto*360 

t  jd i i  r . d 


%  .  1  w  i  r,  i 


is  Drl  mise: r. g  • 
I  s  Dr 2  mssir.q5 


THEN  Dron=Dron+ 36 0 
THEN  Dr to*Dr t c+360 
Test  for  tailwind 
Test  for  tailwind 
Heady i no 


Tai 1 wi nd 

Are  Drl  &  Dr2  missing 
Only  Drl  is  missing 
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1  0  ! 

20  1  COMPUTE  WIND  VELOCITY  USING  RPV  PROGRAM  -  MINI  VERSION 
30  ! 

40  A**“Wind  Speed  =  “ 

50  B*=“Uind  Direction  =  " 

60  D*=  “  m/s“ 

70  E*=“  Degrees" 

80  F*=“UIND  VELOCITY" 

90  INPUT  “Input  L=1  to  compute  V  using  Method  1  or  L=2  for  Method  2",L 
100  IF  L=2  THEN  GOTO  Meth 

110  INPUT  “Input  airspeeds  in  m/s  <X  and  Y>“,X,Y 

120  INPUT  “Input  distances<mi  and  tjmes(s)  < Dx , Tx , Dy , Ty > " , Dx , Tx , Dy , Ty 

130  INPUT  “Input  headings  in  degrees  (Dr-1  and  Dr2  > "  ,  Dr  1  ,  Dr2 

140  PRINT  USING  HeadingjF* 

150  Heading:  IMAGE  15X,13A,2/ 

160  GOTO  Anda 

170  Meth:  INPUT  "Input  airspeeds  in  m/s  <X  and  Y)“,X,Y 

ISO  INPUT  “Input  whether  drift  (I  =  1>  or  correction  <  I  =— 1  >  angle",! 

190  INPUT  "Input  angles  in  degrees  <A  and  6>”,A,B 

200  INPUT  “Input  headings  in  degrees  <  Dr  1  and  Dr  2  ) 11  ,  Dr  1  ,  Dr2 

210  A=I*A 
220  B=  I *B 

2 30  PRINT  USING  HeadingjF* 

240  GOTO  Otra 

250  Anda:  CALL  Comp  1 < X , Y , Dx , Tx , Dy , Ty , Dr  1 , Dr 2, V , D i r  ) 

260  GOTO  Printer 

270  Otra:  CALL  Comp2<  X ,  Y ,  A  ,  B ,  Dr  1  ,  Dr  2 ,  U ,  V  ,  D  i  r  ,  I  > 

280  Printer:  PRINT  USING  T i t 1 e : A* , V , D* , B* , D i r , E* 

290  Title:  IMAGE  5X ,  1 3A , ODD . D , 4A , 5X . 1 7A , DDD  .  D , SA , / 

300  PRINT  USING  Out 
310  Out:  IMAGE  5X,6/ 

320  END 
330  > 

340  ! 

350  !  SUBPROGRAMS 

360  ! 

370  ! 

380  SUB  CompKX,Y,Dx,Tx,Dy,Ty,Dr1,Dr2,V,Dir> 

390  ! 

400  Dir*Q 

410  DEG 

420  L=  0 

430  Xg=Dx/Tx 

440  Yg  =  Dy,’Ty 

450  C=Xg-X 

460  D=Yg-Y 

470  V=<C  2+D'2)" .5 

480  IF  C=0  THEN  C=. 001 

490  IF  D=0  THEN  D=.001 

500  ! 

510  !  computation  of  Direction 

520  ! 

530  G=ATN<ABS<D/C>> 

540  IF  <C>0)  AND  <D>0>  OR  <C<0)  AND  < D< 0  )  THEN  L=1 
550  Dron=Dr1 
560  Drto=Dr2 

570  IF  <  Dron<  =  90  )  AND  (Drto>  =  270>  AND  <Drto<»360  >  THEN  Dr on  =  Dr on  +  36 0 

580  IF  ( Dr  to<  =9  0  )  AND  <Dron>=270>  AND  <Dron<=360>  THEN  Dr t o =Dr to+ 36 0 

590  Dir=Dr1-G  !  Headwind 

600  IF  <  Dr  to<  Dron  >  AND  <L  =  0.>  THEN  Dir=Dr1+Q 

610  IF  < Dr  to>Dron  >  AND  <L«1i  THEN  Dir=Dr1+G 

620  IF  C)0  THEN  Dir=Dir+180  1  Tailwind 

630  IF  Dir  >36  0  THEN  Dir=Dir-360 

640  IF  D l r <  0  THEN  Dir=Dir+360 

650  SUPFND 
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o  6  0  1 

670  ! 

oSO  SUB  Comp2<  X  ,  Y  ,  H  ,  6  ,  C'r  t  ,  Or  2  .  U,  V ,  D  i  r  ,  I  > 

»'*  -ft  o  I 

7uO  D l r  =  0 

7  la  K  =  0 

720  u -•  o 

730  DEG 

.40  c=x*thn<h> 

750  D=Y*THN<6> 

76  0  V=C C " 2+D  *  2  >  * . 5 

770  IF  C=G  THEN  C=.GG1 

730  IF  D=o  THEN  D=.001 

79  0  ! 

000  I  Computation  of  Direction 

310  i 

820  F'  =  H  T  N  <  HB$<  C/O  )  > 

830  IF  <.C>0  >  AND  <D>0>  OR  <  C<  0  >  HND  <  D<  0  )  THEN  K=1 
94  0  Dr on-Drl 
850  Drto=Dr2 

860  IF  <0ron<=90)  HHD  <Drto>=270>  HHD  < Dr t c<=36 0 >  THEN  Dr on-Or on* 36 0 

8.;0  IF  <Drto<=9u>  HND  <  Dr  on  >=2  0  >  HMD  <  Dr  cn< -36  0  :>  THEN  Dr  F  o  =Dr-  tc+  36  0 

980  IF  <0>0  >  HND  <  Dror>Dr  to  >  THEN  1.1-1  >  Tc;t  f...r  tailwind 

890  IF  < D< 0 >  HND  <  Orcn'Drtc  I  THEN  0-1  !  Teat  for  tailwind 

900  Dir “Or  1 -p  !  Headwind 

910  IF  < Dr on >Dr  to)  HND  <  K - 1  )  THEN  Dir=Or1+F 

920  IF  <0ron<Drto>  HND  <  K  -  0  >  THEN  Dir==Ot  1  +F 

9  5  0  IF  IJ=1  THEN  0 1 r  =  D  :  r  +  1 S 0  l  Tailwind 

-40  IF  Dir >360  THEN  Dir=Dir-360 

-50  IF  Dir<0  THEN  Oir-=Dir+36G 

460  SUB END 
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ESTIMATION  OF  CEIL  INN  tSTuREO  hS  "RPVC 


OEG 

D  h  T  A  U,  0 

R  E  a  D  X  .  2 

Ot  = "  Ce  j  1  l  ng  -  " 

Bf  = "  meters" 

INPUT  "Set  1  =  1  it  distance  known,  1=2  it  altitude  known..  1=3  it  both  known" 


INPUT  "Input  upward  ele  at  ion  angle  v  Beta  >  in  degrees" , Beta 
ON  I  GOTO  Ini , In2, In3 

Ini  :  INPUT  "Input  distance  •:  X )  lit  meters'  t  depression  angle  < Alpha!  in  degre 
.  X, Alpha 
GOTO  Compl 

In2'  INPUT  "Input  altitude  <  2  >  in  meters  i  depression  angle  <  Alpha  in  degre 
,  2  .  A 1  p  h  a 
GOTO  Comps 

Ir.3  :  INPUT  "Input  distance  «  X  >  &  altitude  <  2  .»  in  meters"  .  X,  2 
GOTO  Comp 3 

Compl  ;  L=X*‘.  TaN*..  Alpha  >  +  ThN-  Beta  1  > 

GOTO  Out 

Comp  2-  C*2*t  I  ♦TaN«.  Beta  '."fwln  wlpti-t  »  > 

GOTO  Out 

Comp  3  ;  C=2  *X  *TAN<  Beta  . 

Out;  PRINT  USING  T  1 1  le  ;  H*  .  ,  El 

Title:  I  MAGE  &  .  5K ,  1  Oh  .  5f«.  2A,  // 

E 


© 

1  0  ' 


2  0 

t 

ESTIMATION  OF  EXTINCTION  .'OEFFK 

I  ENT  STORED  AS 

“  PAOPF  V"  .• 

3  0 

1 

USING  TWO  SLANT  PATHS 

4  0 

1 

5  0 

1 

INPUT 

6  0 

1 

ro 

DATA 

0,0,0 

30 

P£*C< 

X 1  ,  X2 , 2 

90 

"  T 

INPUT 

"Set  1=1  if  distances  known,  1=2 

if  altitude  known. 

1=7  if  both  known 

\  0  0 

ON  I 

GOTO  In  1 , I n2 , In3 

1  1  0 

Ini  : 

INPUT  "Input  values  of  distances 

■,  XI  .  X2  where  X  1  <  X2  > 

in  meters'.  X  1  X2 

1 .5 0  GOTO  Angle 

i  ">  0  Ini’:  INPUT  "Input  value  of  altitude  f  2  >  in  meters'  .  Z 

'  -1  0  G  0  T  0  Angle 

1  5  0  In3 :  INPUT  "Input  distances  <  XI  .5-2  where  XI  X2  :•  and  altitude  •  2  '•  in  meters" 

.  XI  .  .<2,2 

it>i)  GOTO  In 

'70  Angle:  INPUT  "Input  angles  •  81  .  82  where  81^62  •  between  vertical  and  slant 

paths  in  dear ees " . 8 1 , 82 

1 80  In:  INPUT  "Input  slant  path  radiances  <ftl.Fi  wtiera  ft  l  >R2 ;  in  any  standard  u 

■  its"  .  R  1  .Pi 
(AO  ! 

2 0  0  1  COMPUTAT I  ON 

2  1  0  i 
220  DEG 

13  0  DIM  h*  C  32 ] 

24 o  XI =  X 1 /I  000  1 

25  0  X2-X2  -'1  000  1  Charge  units  to  Hlonsterj 

260  2=2/1000  1 

2/o  A*=" Volume  extinction  coefficient  =  " 

280  B*="  1/km" 

'90  IF  I  >3  I  HEN  GOTO  Alt 
5 u 0  61=ATN< X1/Z  ' 

31  0  82  =  ATN’:  X2/2  > 

320  Alt  :  A  1  ph  a  =  LOG<.  R2/P  1  > 

530  IF  2=0  THEN  GOTO  hist 

340  A  2  ph<a  =  A 1  pha/<  Z*k  I.COStB)  >  —  1 . C  0  9-  •  82  >  >  < 

350  GOTO  Out 

560  Oist:  IF  BIO  0  THEN  GOTO  Set 

570  Z=X2/TAN<:  Bi  > 

38  0  B  1  =  A  TN1'  XI  /Z  ' 

390  Bet:  IF  620  0  THEN  GOTO  Go 
400  2=X1/TAN< BI  > 

410  B2  =  ATrk  X2/2  > 

420  Go:  IF  81  <5  THEN  GOTO  'small 

4  7  0  A  1  p  h  a  =  A  1  p  h  a  /  <  X 1  /  S I N 8 1  1  -  X  2 . '  S I N  <  B  2  1  1 
440  GOTO  Out 

45  0  Small:  A 1  pha  =  A  1  pha.  <  X2*>  1 . 'T  AN1  8.  ■-  1  /S  I  N<  B2  >  .■  1 

46  0  ! 

470  !  OUTPUT 

4  8  0  1 

490  Out:  PRINT  USING  Title: At. Alpha, Bf 

5  00  Title:  IMAGE  ©,5X  32a  ,  DD  .  ODDO  .  nA  .  // 

510  END 
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pfvolt 


© 

1 0  1 

20  '  ESTIMATION  OF  EXTINCTION  COEFFICIENT  -STORED  AS  “ 


3  0 

4  0 

i 

i 

USING 

VOLTAGES  AS  INPUT 

5  0 

6  0 

i 

1 

INPUT 

70 

DaTa 

0,  0 

8  0 

READ 

X  ,  2 

9  0 
.  I 

)  0  0 

INPUT 

"Set 

1  =  1  if  distance  known..  1=2  if  altitude  known.  1  =  3  if  both  known 

ON  I 

GOTO  Ini , In2 , In3 

1  t  0 

Ini  : 

INPUT 

"Input  value  of  distance  i.  X  >  in  meters''./ 

1  2  0 

GOTO 

Angle 

1  30 

In2  : 

INPUT 

"Input  value  of  altitude  > Z >  in  meters", 2 

1  4  0 

GOTO 

Ang  1  e 

150 

I  n3  : 

INPUT 

"Input  distance  <  X  >  and  altitude  • 2  '■  in  meters" ,  X  ,  2 

1  6  0 

GOTO 

I  n 

170  Angle:  INPUT  “Input  angle  CBeta'  between  vertical  and  el  ant  path?  in  deqre 

6s" .Beta 

130  In:  ChLL  Vo  1 1< R 1 , R2 , V 1 , V2 , A , B i 
13  0  1 

000  1  COMPUTATION 

210  1 
220  DEG 
23  0  DIM  A$C32J 

240  X=X/1000  1  Change  units  to  kilometers 

25  0  2  =  2/ 1  0  00  ' 

260  A*  =  "Vo  1 ume  extinction  coefficient  =  ” 

2  70  B*="  f  /km" 

230  IF  1  =  3  THEN  Beta=ATN<  X/2  > 

230  Alpha=L0G- R2/R1  ) 

300  IF  2=0  THEN  GOTO  Diet 

•■1  0  A  lpha=A  1  ph  a  A  Z*<  1  - 1  /COST  Bet  a  -  >  - 

320  GOTO  Out 

330  Oist:  Alpha=Al  pha/(  X*f  1 /TAN- Be  t  a  >- 1 IN<  Be*  a  >  •  • 

74  0  1 

350  '  OUTPUT 

36  0  1 

.3  7  0  Out:  PRINT  USING  T  i  1 1  e  ;  A*  ,  A 1  ph  a  .  B  * 

33  0  Title:  IMAGE  3 , 5X  .  32  A  .  DC'  .  DODO  .  6  A  .  ' 

390  END 
4  0  0  1 

410  1  SUBPROGRAM  TO  CONVERT  VOlThGE  TO  Rwi-i  hmCE 

420  1 

430  SUB  Volte Rt ,R2, Vt ,V2,A,B > 

440  INPUT  11  Input  voltage  coefficients  -h.B-"  .  a.B 

450  INPUT  “Input  voltages  for  vertical  and  slant  paths  •  VI  V2  "here  V ;  VI  . 

•  J  T> 

r  c- 

460  R 1 =A  +  B* V  t 
470  R2=A+B*V2 

480  SUBEND 
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EOT  i  1'Ih  T  i  C'H  OK  £■•  '  i  H>:  r  I  Cn« 
USING  TiV.i  FhIH1:  hh 

h>  IhplIT 


l  H  T  H  'I  lj  ,  0 

F  t  H  0  s)  .•w.r 

I  ,tF  1 1  T  ■  x4  t  i  !  ,  «  J  1  j  • 

■jH  I  -PJ  T 0  lr,  l  :  r. .  :  r-  ~ 

il  l.  i  NFiO  i  '  1  i.-  -i  r  3  .  *  • 

0  0  r  0  ii  '  "3  1 

In*:.  INPUT  "Si. put  ■■.!•- 
ijOTi.'i  Hr.glii 

l •  :  INPUT  "In-  .  t  !  l  J  '  -i 

GOTO  In 

IMP  O’  Ir.p  -.*•  • 

■i  T  in  do  '3'  f  -  3  '  f  ■  t  . 

I  fi  .  ’  Hi.  L  •'  '.If  "  '  K. 

I  E  o  IT  F’U’mT  I  I  j  N 


I.'IM  h  I E  Tv,  1 
l  =  i  ll)  i.i  it 
j  =  i  i.i  0  (i 

■=;  i  ij  n  (.1 

H«=  “Vt  I'JITi--  •:  t  iff-T  .  .  •  .  • 

~  i  -  "  I  i.  Ill  " 

I P  I  V  THEN  N"t0  ill 

p  i  -  h  r  n  ■  >.i.  r  • 

t^^rui  :<o.T  i 
hii  H  I  p  h  s  - 1  i.i  G  |  P'  i-  P  i  1 
IP  0=J  THEN  'jOT'.i  I'iji; 

H  i  pN  3  -  H  1  p  Pi  3  *  1  I  'Pl'i.  i.  1 

1.,  1.1  TO  Out 

[■lit  I P  t- ’  THEN  '..i'll 

o  -  -.o.  r Hie  go 

Ci  =hT N'  •:  i  Z 

b- t  I  r  GO  ■  ‘HEN  j  .*■  - 

J--M  .  'hN  bit 
GO  -h  T f )  '0 

G  :  IP  El  G  THEN  O'  i T  "••-.1 

h  |  P  t,  3  =  H  l  pt>  =t  1  >  I  N  b  1  “  ■  I  ’  ‘ 

GO  10  Out 

Oiii^ll  h  !  pr.  jui  pr.  j  •  0*  i  T  ,-n  t 


PRINT  '■  I  N'j  r  i  t  ! -r  mI  Hip ' 
IMhOE  ij  *0  .T^h  I’C'  for-fi 


'.t.lp  f  P  i  ii.P  hN  7  0  i.  On  -  t  P  ' 

Oi.il:  /..It  PI  P'0  VI  '0  H.L 

i  r  i  P  lj  T  I  r  i  [.  t j  f  i  t  f'lfi  'I'/.it  Kit. 

InFuT  Ir.pu*  v,!*i.)5i  f  '  t'i'h 
P  I  -  H  <  F  .  <  '/'  1 
p,  -  •  ■  »  P  * V 0 
..  •(■  r  hi 
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